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ABSTRACT 

With  constant  firing,  metal  fatigue  produces  cracks  In  a gun  barrel. 


The  useful  life  of  the  barrel  comes  to  an  end  when  a crack  develops  to  a 
critical  sise.  The  theory  of  Fracture  Mechanics  suggests  a formula  for 
crack  sise  growth  rate.  This  formula  can  be  used  to  determine  the  life 
of  a barrel  depending  on  the  initial  and  critical  crack  iizea  ah^  other 


factors.  The  initial  crack  site  turns  out  to  be  a dominant  factor. 
Unfortunately,  accurate  measurements  are  not  generally  available  on  the 

initial  crack  site.  In  this  paper,  s>  prepese  a simple  probability  model 

V £ yt  ^ o \ r J 

for  the  initial  crack  site  and  this,  in  turn,  leads  to  a probability 

A 

distribution  of  the  life  of  the  barrel.  This  last  distribution  is  the 
well  known  exponential  distribution  with  a location  shift.  The  simplicity 
of  this  final  result  is  one  of  the  factors  that  make  the  model  appealing. 
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A PROBABILITY  MODEL  FOR  INITIAL  CRACK  SIZE 
AND  FATIGUE  LIFE  OF  GOlt  BARRELS 

by 

F.  Proschan  and  J.  Sethuraman 
the  Florida  State  University 


1.  Introduction.  Suppose  that  there  is  a crack  of  sire  b in  a (tun 
barrel  after  N rounds  have  been  fired.  The  rate  at  which  the  crack  sire 
will  grow,  ^ , is  a central  topic  of  study  in  Fracture  Mechanics.  The 
following  formula,  taken  from  Davidson,  Throop  and  Felner  (1972),  represents 
a reasonably  simplified  version  of  many  more  complicated  ones  available. 


(1.1) 


db  C(AK)m 
dN  E Sf  KIc 


where 

AK  " Increase  in  stress  intensity 
E ■ elastic  modulus 
Sy-  " tensile  yield  strength 
Kyc  * fracture  toughness 

m ■ a constant  between  2 and  4 
C • a material  constant. 


Throop  (1972)  gives  a simple  expression  for  AK,  which  depends  on  the  sice 
b and  the  shape  of  the  crack,  as  follows: 

(1.2)  AK  - uS/ifb  , 


where  . 

M *+l 

S ■ bore  stress  - p-= — , 

« -1 

p ■ pressure 

«t  • ratio  of  outer  diameter  to  inner  diameter 

) 1.5  if  the  crack  is  a frontal  notch 
1.0  if  the  crack  ll  semi-elliptical 
0.5  if  the  crack  is  semi-circular. 
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(1.3) 
where 

(1.4) 


db  _ .m/2 

5« ■" b • 


1)  - tCo"s".m/2)/(*  Sy  K,cl 


Suppose  that  the  barrel  is  fired  a few  rounds  to  provide  a heat  check 
pattern.  The  size  of  the  largest  crack  produced  at  this  stage  is  called  the 
initial  crack  size  and  denoted  by  b . The  number  of  rounds.  L.  required  for 
this  crack  to  grow  to  the  critical  crack  size  bc  is  the  life  of  the  barrel. 
From  (1.3),  it  is  easy  to  obtain  an  expression  for  L: 


(1.5) 


L • 


2 ..  -(m-2)/2  _ -(*-2)/2 

(m-2)  D *”o  c 

^ [log  b£  - log  bQ]  .If  m - 2. 


] if  m > 2. 


It  is  easy  to  see  from  (1.5)  that  while  the  effect  of  b£  on  L is 
minimal,  the  effect  of  b0  is  significant.  For  Instance,  an  Increase  from 
4 in.  to  5 in.  in  b increases  L only  by  -r~—  , when  m - 4.  *or  the 

C Z\nJ 

same  m,  the  Increase  in  L when  bQ  is  decreased  from  .01  in.  to  .001  in., 
900 

is  — jj— . Also,  the  initial  crack  size  is  difficult  to  measure  since  it  is  so 

small.  It  would  therefore  be  useful  to  model  the  distribution  of  b . ue  do 

o 

this  in  the  next  section  and  show  how  this  leads  to  a simple  distribution 
for  L. 

To  the  best  of  our  knowledge,  there  have  been  only  two  previous  attempts 
made  to  Introduce  probability  models  in  the  problem  of  fatigue  crack  growth. 
Racicot  (1976)  assumed  fixed  values  for  bQ  and  b£,  and  various  probability 
distributions  for  the  severel  fectors  in  !>.  He  then  obtained  empirical  results 


3 


for  Che  distribution  of  L by  simulation  on  a computer.  Hanagud  and 
Uppalurl  (1976)  used  a formula  for  crack  growth  rate  which  is  slightly 
more  complex  than  (1.1)  and  Involves  a quantity  r which  is  the  ratio 
of  the  maximum  stress  intensity  factor  to  the  minimum  stress  intensltv 
factor.  They  assumed  that  both  AK  and  r were  random,  and  obtained 
approximations  to  the  expected  crack  sise  growth  rate.  They  used  these 
results  to  study  Improvements  in  reliability  obtainable  by  examination 
and  repair.  Since  examination  and  repair  are  not  possible  in  a gun 
barrel  operating  in  the  field,  their  results  are  not  applicable  to  our 
problem. 

2.  The  probability  model. 

The  first  few  rounds  of  firing  produce  a heat  check  pattern.  Let 
this  pattern  contain  N cracks  having  sixes  C^,...,C^.  Then  the  initial 
crack  else  bQ  is  given  by 

bQ  • max  {<^,...,0^}  . 

In  our  probability  model,  we  assume  that  C^,  Cj, . . . are  Independent  and 
identically  distributed  according  to  a uniform  distribution  on  fO,Bl. 

We  also  assume  that  * is  an  Independent  variable  governed  by  a Poisson 
distribution  with  parameter  X.  A physical  interpretation  for  B and  X 
is  as  follows:  B represents  the  maximum  possible  initial  crack  sise  and 
X is  a measure  of  the  number  of  cracks  in  the  heat  check  pattern.  Our 
probability  model  depends  on  Just  these  two  easily  Interpretable  quantities. 


N. 
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Thus, 

P(bQ  Sb)  • P(max(Clt...,CN)  i b|w  5t  1) 


l P(max(C.,...,C  ) i b)  P(N  - n) 
o-l 1 n 


[ P(«  - n) 
n-1 


so  that 


(2.1) 

Xb/B  . 

P(bQ  lb)  , 0 1 b 1 I. 

• - 1 

The  distribution  In  (2.1)  Is  our  probability  model  for  the  initial 
crack  else.  Xt  is  physically  motivated  and  simple.  A sample  application 
of  this  model  given  at  the  end  of  this  sactlon  provides  more  support  for  the 
model. 


If  the  distribution  of  D (for  definition  see  (1.4))  or  the  distribution 

of  the  factors  that  sntar  D ara  known,  one  can  combine  that  knowledge 

with  the  probability  distribution  of  b in  (2.1)  to  obtain  the  distribution 

o 

of  the  Ilfs  L of  the  barrel. 

As  a sample  application,  wa  will  obtain  an  approximation  to  the  distri- 
bution of  L.  Assume  that  T>  is  known  and  m ■ 4.  Formula  (1.5)  states 
that 
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Thus,  for  any  x > 0, 

'j1  s(dS  ‘ 5E^)  + xdbJ 

■ p{  r * » + f)'1} 

*x-fl 


(2.2)  * 1 - e’x, 

where  - stands  for  approximately  equal  to  and  ve  have  used  the  relations 
(1  + ~)_1  - 1 - j and  eX  - 1 * eA  , 

which  are  good  approximations  even  for  moderate  values  of  X.  The  conclusion 
in  (2.2)  may  be  summarised  as  follows: 

If  the  number  of  all  crack  slses  in  the  heat  check  is  moder- 
ately large,  the  life  of  a gun  barrel  has  an  exponential 
distribution  with  a location  shift. 

In  the  general  case,  where  h j»  assumed  to  be  random. 
the  distribution  of  the  life  of  the  barrel  becomas  a mixture 
of  shlftsd  exponential  distributions. 

The  assumptions  of  a uniform  distribution  for  initial  crack  also  and  a 
Poisson  distribution  for  the  number  of  cracks  in  the  heat  cheek  pattern  are 
not  required  for  a result  like  (2.2)  and  were  made  here  just  to  Illustrate 
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the  proof.  In  fact,  it  can  be  shown  from  the  theory  of  limiting  distributions 
for  the  maximum  (or  minimum)  of  large  numbers  of  random  variables  that 
the  distribution  of  L can  be  approximated  by  a shifted  ”eibull  distribution. 
In  the  general  case  where  h is  assumed  to  he  random,  the  distribution  of 
L can  be  approximated  by  a mixture  of  shifted  Weihull  distributions  with 
a common  shape  parameter. 

It  is  hoped  that  currently  available  data  on  gun  barrels  will  be 
studied  further  to  test  the  validity  of  our  model. 

3.  Conclusions. 

(1)  A simple  probability  model  for  the  initial  crack  size  is  given  in  (2.1). 

(2)  It  is  derived,  as  a consequence,  that  the  life  of  a barrel  has  a 
shifted  exponential  distribution. 

(3)  Under  practically  no  assumptions,  it  can  be  shown  that  the  life  of 
a barrel  has  a shifted  Weibull  distribution. 

(4)  When  D as  defined  in  (1.4)  is  assumed  to  be  random,  the  distribution 
of  the  life  of  a barrel  is  a mixture  of  shifted  Weibull  distributions  with 
a common  shape  parameter. 
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